Cavity driven by a single photon: conditional dynamics and non-linear phase shift 
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We apply the stochastic master equations (quantum filter) derived by Gough et al. 1-3 to a 
system consisting of a cavity driven by a multimode single photon field. In particular, we analyse 
the conditional dynamics for the problem of cross phase modulation in a doubly resonant cavity. 
Through the exact integration of the stochastic equations, our results reveal features of the problem 
unavailable from previous models. 
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I. INTRODUCTION 

Single photons constitute an important element for a 
variety of quantum technological applications. Some- 
times referred as a "flying qubit" , propagating single pho- 
ton wavepackets can serve, for example, as the qubit s 
used for linear-optical quantum computation [4 , or as 
carriers of quantum information between nodes in a 
quantum network [5H7J. 

These possibilities have triggered an enormous exper- 
imental effort toward the production of single photons 
on demand [8]. Theoretically, the focus has been on the 
development of a proper description of the continuous- 
time multimode single photon field and its interaction 
with other quantum systems, such as atoms or cavities. 
In Ref. [9], for example, Gheri et al. introduced a for- 
malism to describe the dynamics of a quantum optical 
system driven by a single photon field in terms of a gen- 
eralised master equation. A proper treatment for the 
multimode single photon field has also been taken into 
account in some photonic quantum gates [TO] [TT] and 
also in the analysis of the loading probability in trapped- 
atom quantum memories [12] . 

More recently, Gough et al. [TH3] developed a stochas- 
tic master equation (SME) [T3l[T4] , or quantum filter [T5l - 
[T7] to describe the situation where a quantum system is 
continuously monitored and driven by a variety of non- 
classical states, including single photon fields. This de- 
scription is important not only to analyse the conditional 
dynamics of systems undergoing continuous measure- 
ment but is also an essential step to model measurement- 
based feedback. 

Quantum feedback control has been recently used in a 
range of problems from the preparation and protection 
of quantum states [T8U20] to the stabilisation of Bose 
Einstein condensates [2TJ [22] . In the same way that feed- 
back and other control techniques became ubiquitous in 
classical technological applications, we expect that they 
will play a similar role in future quantum enabled tech- 
nologies. In a network of quantum systems, for example, 
where one could imagine that some of the channels con- 



necting nodes will be fed back to the system while oth- 
ers will be measured, it becomes essential to provide a 
physical description of the conditional dynamics of such 
processes. 

In this paper we apply the formalism developed in pj[3] 
to describe the conditional dynamics of the cross phase 
modulation in a doubly resonant cavity. In this prob- 
lem, a single photon pulse interacts with a coherent field 
through a cavity nonlinearity inducing a phase shift in 
the coherent mode. Our goal here is to investigate how a 
proper stochastic model and its solution can bring new in- 
formation about the conditional phase shift process. The 
paper is organised as follows: In Section[ll]we present the 
stochastic Schrodinger equations (SSE) corresponding to 
photodetection and homodyne measurement of the single 
photon channel. This is an alternative to the SME de- 
rived in [T]-[3] and an explicit construction of the mas ter 
equation unravelling proposed in [9 . In Section 



III 



we 



analyse the simple model of a driven single mode cav- 
ity. The simplicity of the system allows us to develop 
intu ition about the driving and the monitoring. In Sec- 
tion IV we discuss the model [TTj of a doubly resonant 



cavity that is driven by a single photon pulse and a coher- 
ent field. We analyse the problem of the non-linear cross 
phase shift imprinted by the single photon on the coher- 
ent field conditioned on the measurement of the single 
photon mode. Finally, in Sec. [V] we conclude and discuss 
future perspectives. 



II. THE MODEL: CONDITIONAL EVOLUTION 

We start this section with a brief summary of some of 
the results from p~H3] to establish the notation and pre- 
pare the ground to present the SSE for a single photon 
channel. The general scenario we want to describe is the 
one shown in Fig. [T] where a continuously monitored sys- 
tem is probed by a single photon field. The state of this 
multimode field is given by |1^> = / °° dt£{tp(t)\0) [23], 
where / °° dt|f(t)| 2 = 1. 
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Before moving into the monitored evolution, it is in- 
structive to look at the unconditional dynamics given by 
the master equation. Let p be the density matrix for the 
joint state of the field and the system. Now, in order 
to describe the dynamics of the system itself one needs 
to trace over the field variables. When the field is in 
a single photon state, the relevant unconditional system 
density operator is given by pu(t) = (l^\p\l^). The mas- 
ter equation for the system derived in [THSJ [5] is given by 



gle photon state is that the dynamical equation for the 
system operator pn(t) depends also on the operators 
Poo ? Poii an d Pio- This is in contrast with the situa- 
tion where the field is in a Gaussian state and there is 
only one equation for the system operator p s = Trfi e id [p] . 
Note, however, that pn(t) is the only physically rele- 
vant operator since the expectation value of any sys- 
tem operator X s can be calculated directly from it using 
(X s )(t)=Tr[X sPll (t)}. 
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where the Liouvillian term Cp(t) — —i[H, p] + V[L]p con- 
tains the Hamiltonian and the decoherence part in the 
usual Lindblad form V[L)p = Lptf - l/2(tf Lp + ptf L). 
An important consequence of the field being in a sin- 
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detection 
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FIG. 1: Schematics of a quantum system probed by a single 
photon field. The dynamics of the system will be conditioned 
on the measurement outcomes Y(t). 

Now we want to describe the dynamics of the moni- 
tored system. From the general filtering equations de- 
rived in [T]-[3] , we can write the conditional master equa- 
tion for the system. In the case of homodyne monitoring 
they are 



d Pll (t) = [£ Pll (t) + \poi(t),lS]£(t) + [L,pio(t)]rW] dt 





+ [Lp n (t) + 


Pii(f)£ f + Pw(t)C(t) + fln (*)£(*) - K t pu(t)] dW(t), 


(2a) 


dpoi (t) 


= [Cpoi(t) + 


[L,p 00 (t)}C(t)} dt + [Lp 01 (t) + p 01 (t)tf + p 00 (t)C(t) ~ K t p 01 (t)] dW(t), 


(2b) 


dp 10 (t) 


= [Cpi (t) -f 


- [poo, £*]£(*)] dt+ [L Pl0 (t) + Pl0 (t)tf + p 00 (t)Z(t) - K tPl0 (t)] dW(t), 


(2c) 


dp 00 (t) 


= £p 00 (t)dt 


+ [L Poo (t) + p 00 (t)L ] - K tPoo (t)] dW(t), 


(2d) 



I 

with dW(t) a Wiener process obeying dW(t) = dY(t) — For the case of a photocounting monitoring the equations 
K t dt and K t = Tr [(L + L^)pn + poi£(i) + PioC(t)] • become [3] 

I 



dpn(t) 


= [Cpu(t) + 
+ K l (Lpn 


[p 01 (t),L^(t) + [L,Pio(tW(t)]dt 
(t)Lt + Lp w (t)C(t) + Poi(t)L^(t) 4 


-p 00 \m\ 2 )-pu(t)}dN(t), 


(3a) 


dpoi (t) 


= [Cp 01 (t) + 


[L,p 00 (tW(t)]dt+ [vt 1 (Lp 01 (t)tf 


+ Lp 00 (t)Z*(t))-p 01 (t)] dN(t), 


(3b) 


dpw(t) 


= [Cp 10 (t) + 


[poo, £*]£(*)] dt + [vt 1 (L Pl0 (t)tf + 


p 00 (t)L^(t)) ~ PioW] dN(t), 


(3c) 


dp 00 (t) 


= Cp 00 (t)dt- 


f [v^LpooftlS - pooW] dN(t), 




(3d) 
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where dN(t) is a compensated Poisson process obeying 
dN(t) = dY - v t dt, with v t = Tr[tf Lp n + LpioCCO + 
L^poi£(i) + Poo IC(^) 1 2 ] - Note that by taking the classi- 
cal ensemble averages over the noisy processes we obtain 
back the master equations Eqs.Q. 

A crucial aspect of the single photon master equations 
is their non-Markovian behaviour. This is important as 
SSE unravellings of the master equation are only guar- 
anteed to exist if the later is Markovian. The solution 
to find the SSE corresponding to Eqs.Q and (|3| is to 
embed the system in a larger Hilbert space that contains 
an ancilla which generates the input field state. Such 
Markovian embedding has been suggested in [9] for the 
case of a single photon field and derived in [2j [3] for other 
non-classical input fields in the context of quantum fil- 
tering theory. 

Finding this ancilla for general input states is a hard 
task but can be simple for a single photon field. Gheri 
et al. [9 proposed a cavity initially with one photon as 
the auxiliary emitter while Gough et al. [2j |3] suggested a 
two-level system initially in the excited state. Using the 
latter as the choice for our ancilla and using a cascade 
approach we can write a master equation for the total 
system (cavity + ancilla) as 



where 



d j- t =-i [H t ,q]+V[Lt] g, 



and 



H'-p — H 



(r(f)<r+a-£(t)aV_), 



(4) 



(5) 



(6) 



where a + = (cr_)t = \e){g\ and \e)(\g)) represent the ex- 
cited(ground) state of the two-level ancilla. As already 
suggested in [9 j, from here it is easy to write quantum 
trajectory unravellings as Eq.Q is Markovian and in the 
Lindblad form. Explicitly, the stochastic Schrodinger 
equations for homodyne and photocounting are, respec- 
tively 



d|*> = \-iH T +{L ] T )L T - 
+ (L T -(L T ))\*)dW, 



LlpL' 



t ±j T 



(L T )(L T ) 




and 



d|*> 



Li 



T lJ T 



|#)dt 



(l1j.Lt) 



1 \V)dN, 



(8) 



with the point process dN obeying dN 2 = dN and 
E[dN(t)] = (L ] T L T )dt. The connection between Eqs.{7|) 



and ([2| and Eqs.Q and ([8| can be established using the 
correspondences [2, 3||9] 



Poo 
Poi 
pio 

pll 



w(ty 

Qeg 

Qge 

Qee H~ Qggi 



(9) 
(10) 

(11) 
(12) 



with w(t) = J t °° \£(s)\ 2 ds and 6ij = 

The advantage of the SMEs Q and ([3| is that they 
describe the conditional dynamics in terms of system 
operators only. The price paid is the necessity of extra 
coupled equations for the auxiliary quantities poo> Poi 
and pio. On the other hand, the description in terms of 
pure states is given by a single stochastic equation 
or (J8|, depending on the choice of measurement process), 
but requires extending the system with the addition 
of an auxiliary 2-level system. Numerically, the SSEs 
have the advantage of simulating at the state vector 
level requiring 2 x TV levels (with TV the dimension of 
the system of interest) instead of the 4 x N 2 for the SMEs. 



III. SIMPLE EXAMPLE: CAVITY DRIVEN BY 
A SINGLE PHOTON 

We will now consider the situation where the system 
of interest is a single mode cavity with decay rate n. In 
this case, the conditional dynamics of the driven cavity 
is given by the models described in the previous section 
with L = yfHa, where a is the annihilation operator for 
the cavity field. To start with, let's analyse the uncon- 
ditional evolution and look at the dynamics of the mean 
number of photons in the cavity assuming that it starts 
in the vacuum. Averaging over the noise processes in 
Eq. ([2| or ([3| we can write the system of equations nec- 
essary to describe the evolution of (n)n = Ti[pna^a] as 



d(n)u 

dt 
d(a) oi 

dt 



-K(n) n - V£«a) ir (*) - <a)oiCW)](^) 



-7j( a )oi 



(14) 



where for an arbitrary system operator c we define 
(c)ij = Tr[(pij)tc] and we have at the initial time 
(a)ij(0) = (n)ij(0) = 0. For the single photon shape 
we choose {(t) = y/je-^-^^Q^ - t ), with G(t - t ) 
being the Heaviside step function, which corresponds to 
a single photon emitted by a two level atom with de- 
cay rate 7. For this case, the exact solution nu(t) = 

ArfKe^-V (e^-^-^ - G(t-t )/(~f-K) 2 is shown 

in Fig{2]-a for different values of The larger the ra- 
tio 7/^, the faster the photon enters the cavity but the 
maximum number of photons in the cavity is obtained 
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when the single photon and cavity rates coincide, i.e. for 
j/k = 1. Note that Eqs.(13) and (14) are identical to 



the ones corresponding to the situation of a cavity driven 
by a coherent field with amplitude e(t) = i^/n^{t). This 
means that just by looking at the average number of pho- 
tons in the cavity one cannot distinguish between a single 
photon or an equivalent coherent driving. 

Obviously this does not hold true for higher order mo- 
ments as the distributions for the field inside the cav- 
ity are different in each case: for the coherent driving 
the vacuum is displaced but the state remains coherent 
while the single photon driving mixes a single photon 
state with the vacuum component of the field as shown 
by the Wigner function plot in Figure [3J 



a) 



0.6 r 




FIG. 2: Average number of photons in the cavity for a sin- 
gle photon driving, a) Exact master equation solutions for 
j/k = 1 (blue, solid line), 10 (black, dot-dashed), and 0.1 
(red, dotted), b) Single trajectories for photodetection (top, 
red curve) and homodyne (solid black line) monitoring with 
j/k = 1. The exact master equation solution is also shown 
(blue dashed line) for comparison. 



This difference becomes evident when the system is 
being monitored. For the photocounting monitoring the 




FIG. 3: Wigner function for the field inside the cavity cor- 
responding to the single photodetection trajectory shown in 
Fig. [2}b at nt = 2.8. Just before the jump the field is in a 
mixture between the vacuum and a single photon state. 



set of equations needed to describe the system is 

d(n) n = [-K,(n)u - V^(a) if*(t) - >/*>(<>>} oi€(t)] dt 

+ [^t -1 WooieWI 2 - Mil] dN(t), (15a) 
d(a) 01 = ~(a) i - y/K(l)oo€(t) dt 

+ v^(rc> 00 £(*) - <a>oi] dN(t), (15b) 
d(l)oo = [^ -1 *»oo - Woo] dN(t), (15c) 
d(n) o = -K(n)oodt - (n) 00 dN(t). (15d) 

Note that, because we assume that the cavity starts in 
the vacuum, at any given time there can be at most one 
photon in the cavity and we neglected expectation val- 
ues with higher order in a and a) such as (a)a)a)ij in 



Eqs.(15). Note also that (n)oo(0) = and only the first 
three equations need to be considered to examine the 
conditional number of photons in the cavity. For the co- 
herent driving the equation is 

d(n) = (-tt(ra) +i(e(t)(a) - e*(a f )) dt 

+ (-(n) + (a t2 a 2 )/(n))^7V. (16) 

In this case, the stochastic term cancels (the state is an 
eigenstate of the jump operator a), showing that the cav- 
ity field is not affected by detection events. The evolu- 
tion therefore coincides with the unmonitored one and 
is shown by the dashed line in Fig(2]-b. The conditional 
dynamics for a single photon driving is quite different: In 
the case of a homodyne monitoring the single photon en- 
tering the cavity is perturbed by the measurement noise 
and ends up evolving back to the vacuum state, while for 
photodetection a detection induces a jump back to the 
vacuum state (see Figj2]-b). 
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IV. TWO-MODE CAVITY AND NONLINEAR 
PHASE SHIFTS 



The results from the previous section will help us anal- 
yse the situation shown in Fig. [4j A doubly resonant 
cavity containing a Kerr nonlinear medium is driven by 
a single photon in one mode (mode a) and a coherent 
field in the other (mode b). Due to the nonlinear inter- 
action between the modes, the presence of a photon in 
mode a will induce a phase shift in mode b. This nonlin- 
ear cross-phase modulation is the basis for the proposal 
of all-optical two-qubit gates [24] and has been recently 
analysed from the point of view of multi-mode propagat- 
ing single photon fields [TO] [TTJ [25]. The Hamiltonian 
describing the system is 



H = x^ba^a + S a 



(17) 



with x being the nonlinearity, S a the detuning of mode 
a and f3 the coherent field driving mode b. Here we will 
examine the dynamics of this system when only the mode 
driven by the single photon field is monitored. While this 
will shed some light into the conditional dynamics of the 
system, our aim here is not to explore in depth the cross 
phase modulation problem but rather illustrate the use 
of the quantum filtering equations described in Sec.[ll|in 
a more complex scenario. 



cavity 
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FIG. 4: Schematic representation of a doubly resonant cavity 
containing a Kerr nonlinear medium. The cavity is driven by 
a coherent field in one mode and by a single photon in the 
other. The output of the single photon mode is monitored by 
a photodetector while the other mode remains unobserved. 



To describe this system we use the stochastic equa- 
tions presented in Section [TTJ Note that because only 
mode a is being monitored, operators for the b mode 
only appear in the Hamiltonian and in the decoherence 
term. We then simply need to replace L = ^Jn^a and 
Cp = -i [H, p] + V\^fK~ a a\p + V\^fk~ h b\p in Eq. Q to ob- 
tain the exact conditional dynamics under photodetec- 
tion. The equation for pu(t) (|3^i) can be written explic- 



itly as 
dpn(t) 



-i [H, p] + V[^T a a]p + V[^b]p] dt 
[V^lpoi (t), a^(t) + y/iQa, p 10 (*)] dt 
v~ x ( ft a apii(t)a f + V^ a Pio(£)<T (t) 
/^Poi(t)a^(t) + poo\m\ 2 ) - Pn(t) ] dN(t). 

(18) 



The changes in mode b due to the single photon in mode a 
can then be obtained by writing the equations of motion 
for the expectation values of the field b. Using Eq. (18) 



(and the equivalent ones for the other pij components) 
together with the relation d(b)u = Tr[bdpu], we get 



d(b)n = -ifi - ix(bn a )i 



Kb(b)i 



dt- 



((bn a ) n + (ba) 10 C + + Woo|£| 2 )] dN. 

(19) 

Similarly to the single mode case, the equation for (6)n 
couples to other moments and the system of equations 
can get cumbersome specially in the conditional case. For 
this reason we omit these equations and will focus now 
on the numerical simulations of the stochastic equations. 
The introduction of the extra mode b as compared to 



the situation analysed in Sec. Ill increases the dimen- 
sionality of the problem and suggests that it could be 
more efficient to use the SSE (|| instead of the SME 
However, the model that we propose introduces an ex- 
tra difficulty which lies in the fact that mode b is not 
being monitored and therefore the system needs to be 
described by a density matrix. In this case it is still pos- 
sible to use a description in terms of pure states via Eq. ([8| 
but it is necessary evolve an ensemble of pure states with 
appropriate weights [26]. For the parameters we used, 
simulations using the full SME turned out to be more 
efficient due to the excessive number of paths needed for 
convergence using the SSE. 

For our simulations we fixed the decay rate of mode 
a and set 7 = n a so that we have the best matching 
between the cavity and the incoming single photon as 
explained in Sec. |III[ We also assume that, as the sys- 
tem is monitored, we estimate the number of photons in 
mode b and feedback this information to adjust mode a 
detuning as S a (t) — —x( n b(t)). This ensures that mode 
a is always resonant even with the nonlinearity present, 
as can be seen from the unconditional equation for (a)n 



d(a) 



11 



dt 



-iS a (a)ii-ix(rib a) 



K a (a) 



11 



11 o V^fWWoi. 

(20) 

The phase shift imprinted in mode b is proportional to the 
product of the nonlinearity x an d the amplitude of the 
coherent field /?. For realistic cavity parameters the non- 
linearity is small so a large (3 is needed for the change in 
the b field to be appreciable. To be able to see the effects 
and yet have simple simulations, we fixed the nonlinear- 
ity at a relatively large value x = ^o/10 so that the basis 
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describing mode b could be truncated at a numerically 
reasonable level. 

Fig. [5] shows the X5 = Re((6)n) quadrature of mode 
b as a function of time for different values of ^5. This 
quadrature is representative of the change in mode b as 
the orthogonal quadrature P5 doesn't change appreciably 
during the process. For each value of we adjusted 
the coherent driving to keep constant the ratio A/3 = 
4/3x/^5, which corresponds to the maximum conditional 
change in mode b found in [TT] for Kb ^> X- I n their 
model, Munro et al. calculate the phase shift assuming 
that at some random time there will be a single photon 
in mode a (n a = 1) and that it will remain there for 
some random interval T, after which the photon leaves 
the cavity. Because we can't possibly know when the 
photon enters mode a, in our description we only use 
the information about the single photon shape and the 
detection clicks at the cavity output. 




-0.1 2 L 
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FIG. 5: Evolution of the X quadrature of mode b for a sin- 
gle trajectory with Kb/^a = 10 (solid, black), 3 (dashed, 
blue) and 1 (dotted, red). The coherent driving is set to 
P = «f/(4« a ) such that A/3 = 0.1. 

We found that by increasing ^5 not only the time taken 
for mode b to relax back to steady state after a detec- 
tion event decreases, but also the maximum phase shift 
increases. This is in qualitative agreement with the re- 
sults in [TT] which shows that the maximum phase shift 
is obtained for ^ ^> K a . However, the conditional phase 
shifts obtained here can exceed the maximum possible 
value, A/3, found in [TT] . 

To explore how our exact stochastic simulations can 
bring new information about the conditional phase shift 
process, we calculated the maximum conditional phase 
shift for a set of 5000 trajectories and constructed the 
histogram depicted in Fig. [6] for (a) ^ = K a j2 and (b) 
Kb = An a . As the ratio Kb/K a gets larger, the change in 
the b field gets sharper (see Fig. [5| and this is mapped to 
a distribution that favours larger phase shifts as shown 
in Fig. [6}b. An interesting aspect that our model can 
capture is the structure of the distributions in Fig. [6] 
They are divided in 2 parts: one corresponding to small 



phase shifts and therefore to detection events that oc- 
curred shortly after the photon in mode a enters the 
cavity, and another one with larger shifts representing 
detections after mode a had enough time to build up in 
the cavity. The two parts are divided by a gap region 
with values that do not show up in any trajectory. This 
region can be understood by looking at the probability 
of detecting a photon between t and t + dt, which is given 
by Pj — v t dt with 

vt = «a(n)n+v^((a)oirw + («>Siew)+(i>ooiewi 2 - 

(21) 

For a cavity that is not driven, the rate of detection is 
given by the first term in Eq.(|2T|), i.e. it is simply propor- 
tional to the decay rate and the mean number of photons 
in the cavity. When the driving is on, all the terms have 
to be considered to account for the interference that leads 
to the dip in v t shown in the inset of Fig. [6j At that time 
no photons are expected to come out of the cavity as v t 
goes to zero. Note that this behaviour is independent of 
the coupling between the two modes and appears also in 
the single mode problem discussed in Section III A sim- 



ple calculation of the field output intensity also shows 
that this effect occurs for a coherent driving and does 
not rely on the single photon character of the field. 



V. CONCLUSIONS 

We have presented the stochastic Schodinger equations 
corresponding to the quantum filter, or stochastic mas- 
ter equations, for systems driven by single photons. We 
applied this trajectory formalism to investigate the con- 
ditional phase shift of a field that interacts nonlinearly 
with a single photon pulse. By presenting a complete 
quantum stochastic model and numerically solving it, we 
corroborated some of the results in Ref. [11 but were 
also able to predict different effects on the distribution 
of conditional phase shifts. The formalism applied here 
is well suited to model measurement and feedback and 
we expect that the approach will be useful to describe 
other systems such as quantum logical operations, quan- 
tum repeaters, or the interaction of nodes in a quantum 
network connected by single photons. 
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FIG. 6: Histogram of the maximum conditional change in the 
b mode for 5000 different trajectories and (a) Kb = ^a/2 and 
(b) Kb = 4ft a - The inset shows v t for a particular trajectory 
where no photons have been observed up to the time shown 
in the plot. The gaps in the distributions exists because there 
are no photons coming out of the cavity when v t — 0. 
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